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Abstract 

Large eddy simulations of three convective 
turbulence events are investigated and compared with 
observational data. Two events were characterized with 
severe turbulence and the other with moderate 
turbulence. Two of the events occurred during NASA ’s 
turbulence flight experiments during the spring of 2002, 
and the third was an event identified by the Flight 
Operational Quality Assurance (FOQA) Program. Each 
event was associated with developing or ongoing 
convection and was characterized by regions of low to 
moderate radar reflectivity. Model comparisons with 
observations are favorable. The data sets from these 
simulations can be used to test turbulence detection 
sensors. 

Introduction 

Commercial aircraft encounters with unexpected 
turbulence may lead to injuries and fatalities, especially 
among unsecured passengers and crew members. 
Therefore, the National Aeronautics and Space 
Administration (NASA), through its Aviation Safety 
Program, is addressing turbulence hazards through 
research, flight experiments, and data analysis. Recently, 
an initial flight campaign, that used NASA-Langley’s B- 
757, has tested a Turbulence Prediction And Warning 
System (TPAWS). 1 ’ 2 ’ 3 ’ 4 

Aircraft encounters with turbulence may occur 
within either clear air or clouds. Turbulence associated 
with convective clouds, otherwise known as Convectively 
Induced Turbulence (CIT), is a leading hazard to aviation 
operations. Approximately 82% of turbulence encounters 
that lead to injuries are in proximity to atmospheric 
convection. 5 Although aircraft attempt to avoid flying 
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directly into the path of thunderstorms, unsuspecting 
encounters with CIT continue to occur. These 
encounters may occur as aircraft: 1) try to skirt around 
the high reflectivity regions of storms in order to 
minimize the deviation from their planned route, 2) 
encounter unexpected convection that appears to be 
invisible or benign from the aircraft’s radar, 3) encounter 
storm tops that rapidly rise into the aircraft’s flight path, 
and 4) are inadvertently vectored into convection by Air 
Traffic Control (ATC). In some cases of CIT, the ship’s 
conventional radar, which is the pilot’s only onboard tool 
for severe weather avoidance, may not detect any radar 
reflectivity signature. TPAWS is designed to provide 
enhanced CIT avoidance, even in these regions of 
minimal radar reflectivity. 

A second flight campaign to demonstrate the 
feasibility of TPAWS to detect CIT was conducted in the 
spring of 2002. TPAWS is centered on the development 
of radar algorithms that produce a reliable turbulence 
alert of an acceptable time prior to a hazardous 
encounter. These algorithms process standard airborne 
radar data to make predictions, thus requiring no 
hardware modifications to aircraft. Analysis of the radar 
performance during the second campaign is promising. 
Of 49 encounters with moderate to severe turbulence, 
81% of the events were predicted 30 seconds or more 
prior to encounter. 6 

Turbulence detection software must be certified by 
the FAA before being installed into the commercial fleet. 
A major goal of the TPAWS project element is to support 
this certification process. This includes developing a ‘tool 
set’ consisting of atmospheric turbulence data sets, 
hazard estimation models, aircraft performance models, 
and radar simulation software, which used together will 
establish statistically viable radar performance. Hence, 
NASA’s role will be to provide tools and procedures to 
assist in the evaluation of an operational system that will 
result in minimal missed detections and false alerts. 

Presented in this paper will be data sets produced 
from large eddy simulations of three CIT encounter 
environments. Also examined are a hazard estimation 
model and a simple aircraft performance model. The 
hazard estimation model enables predictions of 
turbulence loads based on the vertical wind field, and is 
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Figure 1 . Comparison of the peak load (An) to the computed 5 second variance of the peak load (o An ) for 102 
turbulence encounter cases (Bowles) 6 . 


flight path independent. The output from this product can 
determine regions where turbulence would be a hazard to 
aircraft operations, and will be useful to compare with 
radar simulations that utilized the convective turbulence 
data sets. The aircraft performance model is an 
additional tool that can be used to compare with the 
hazard estimation and radar simulated predicted loads. 
Both models are verified with NASA flight in situ data 
and applied to the simulation data sets. Also, ground- 
based Weather Surveillance Radar - 1988 Doppler 
(NEXRAD) 7 radar data is examined in order to support 
the verification of the simulated data sets. 

Definitions 

Quantification of the turbulence hazard must be 
related to the aircraft response, and should have 
properties that can be estimated by onboard detection 
sensors. Metrics for defining the turbulence hazard can 
be quantified in terms of RMS normal loads (a An ), where 
a An is based on moving 5-second intervals of aircraft 
normal loads (An). 1,2 Using this metric, moderate 
turbulence is defined for: 0.20 g < a An < 0.30 g, and 
severe when a An > 0.30 g. Bowles 6 has shown that the 


peak aircraft normal loads, are roughly proportional to 
a An , as shown in Figure 1. In no event did lAnl ever 
exceed 1 g when a An < 0.30 g. On the other hand, peak 
loads may exceed 1 g in the severe turbulence category; 
thus, turbulence encounter intensities with a An > 0.30 g 
may cause passengers and obstacles in the cabin to 
become weightless, an obviously dangerous situation. 
Therefore, a ‘Must Alert’ criteria has been proposed 
when the operational turbulence detection system 
predicts an event with a An > 0.30 g. Also a ‘May Alert’ 
region is proposed for 0.20 g < a An < 0.30 g, which would 
provide a buffer that would ensure the minimizing of 
false alerts. 

Turbulence Data Sets 

Characterization of the turbulence associated with 
CIT events is critical to the development of new 
technology that would warn for such events. 
Furthermore, these characterizations are critical when 
producing data sets required for the ‘tool set’ that will be 
applied in the certification process. Therefore, three 
particular events of interest, two events from NASA’s 
2002 Flight Tests (Events: 235_05 and 232_10), along 
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with a FOQA* case of severe turbulence encounter that 
occurred near Wilmington, DE are investigated. 
Presented here are discussions of each encounter and the 
corresponding results from large eddy simulations (LES). 

Numerical Model and Initial Conditions 

The numerical model employed for this study is the 
Terminal Area Simulation System (TASS). TASS is an 
LES model with a meteorological framework, and has 
been used in previous studies to examine a wide variety 
of weather/fluid mechanics problems such as microburst, 
wind shear, wake vortices, and turbulence. A description 
of TASS is located in Proctor, Hamilton, and Bowles. 2 

In order to simulate the convective turbulence 
environments encountered during each of the three 
events, appropriate initial conditions must be chosen. In 
TASS, each convective event is generated via 
initialization with the wind velocity, temperature, and 
humidity that represents the ambient conditions of that 
case. Turbulence, at scales important to aircraft response, 
is generated by the dynamic interaction of the convective 
cells and ambient winds. Localized regions of turbulence 
are produced by TASS with realistic intensities and 
accompanying radar reflectivity. Output from the model 
includes three-dimensional data sets of wind fields, radar 
reflectivity, moisture, and temperature. 

NASA Event 235-05 

Observation 

The first event to be presented occurred during the 
second NASA TPAWS flight campaign, which was 
conducted during April and May of 2002. Details of the 
aircraft sensors, flight preparations, and flight 
requirements/restrictions are available in Hamilton and 
Proctor. 1 

This case occurred on May 10, 2002, when the 
NASA B757 encountered turbulence over Charleston, SC 
(Figure 2). During event 235_05, moderate turbulence 
was encountered in nearly cloudless air just above a 
rapidly developing thunderstorm. Radar reflectivity was 
up to 44 dBZ near the flight path but detectable values 
were never penetrated. In fact, the aircraft maintained 
visual meteorological conditions throughout the event, 
skimming what appeared to be light cloud debris for 
about a 2 second period. Along the flight path, no radar 
reflectivity signatures appeared on the ship’s 


Flight Operational Quality Assurance (FOQA) Programs 
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unsafe but otherwise uneventful airline flights in order to 
correct potential problems before they lead to accidents. 


convectional radar, which is purportedly the case in many 
CIT encounter accidents. A data set from the numerical 
simulation of this event could be useful in evaluating 
TPAWS alerting rules, since turbulence may not be 
detected along the flight path due to weak radar 
reflectivity, but may be seen nearby where higher levels 
of radar reflectivity exist. 

The convection that produced the turbulence 
environment was developing along a westward 
propagating sea-breeze front oriented parallel to the 
Atlantic coastline. A number of cells had developed 
along the sea breeze with storm tops between 6 and 8 km. 
The cells were pulsing along the front and dissipating as 
they drifted slowly to the southeast over the relatively 
cooler sea-breeze airmass. Horizontal winds at flight 
level (6.5 km Above Ground Level (AGL)) were weakly 
sheared and storm tops ranged between 6 and 8 km. 

The NASA B757 maintained visual meteorological 
conditions (VMC) throughout the event. At 19:20:39 
UTC, the aircraft passed a relatively high reflectivity (44 
dBZ) region approximately 2 km off the aircraft’s 
starboard side. Nine seconds later the aircraft 
experienced a peak updraft of 1 1 m s' 1 in relatively clear 
air, although some cloud debris surrounded the aircraft 
(Figure 3). Airborne radar suggested radar reflectivity 
levels to be between 16-20 dBZ about 120 m below the 



\ \ \ 

Figure 2. Composite f radar reflectivity (dBZ) and flight 
path with in situ RMS normal loads (g’s) along path for 
Event 235_05. 


1 Composite radar reflectivity is a display of the 

maximum reflectivity from any elevation angle at every 
range from the radar. This product is used to reveal the 
maximum reflectivity above any location. 


3 

American Institute of Aeronautics and Astronautics 




flight path but no apparent reflectivity at the flight path 
altitude. The peak normal loads (An) associated with the 
updraft region were +0.43 and -0.38 g and a An peaked at 
0.22 g at approximately 19:20:50 UTC (Figure 4). The 
aircraft continued to experience successive positive and 
negative normal loads (An) due to a series of updrafts and 
downdrafts in the clear air regions at the lee of the cloud. 

Numerical Simulation 

The simulation for this case included a domain size 
of 20 km x 20 km x 12 km. The total number of grid 
points in the domain is 206 x 206 x 130, with a uniform 
grid size of 100 m above a stretched boundary layer 
region. The initial environmental conditions were taken 
from the Charleston rawinsonde observation at 0000 
UTC on 1 1 May, 2002. The boundary layer for the input 
conditions was modified using Automated Surface 
Observation System (ASOS) data taken from Charleston 
one hour prior to the encounter. 

A brief comparison of the simulated and observed 
convective characteristics is shown in Table 1. A three- 
dimensional perspective of the convective cell is shown 
in Figure 5. Figure 6 shows the simulated radar 
reflectivity at 7.4 km , which compares well with radar 
signatures observed by the Charleston, SC NEXRAD 
radar shown in Figure 7. An altitude of 7.4 km was 
chosen because it had better characteristics with the 
observed event. This is about 900 m above the actual 
flight path. 


In Situ of Event 235 05 



Figure 3. In situ vertical velocity (m s' 1 ) and o An (g’s) 
along path for Event 235_05. 


In Situ of Event 235 05 



Table 1: Comparison of simulated and observed 

characteristics for Event 235_05. 


Variable 

TASS 

Observed 

Peak Storm Tops 

8 km 

6-8 km 

Peak Radar 
Reflectivity near 
Flight Level 

45 dBZ 

44 dBZ 

Cell Motion (toward) 

ESE at 
8.7 m s' 1 

ESE at 
5.3 m s' 1 

Diameter of 
Convective Cell 
(based on 20 dBZ) 

3.5 km 
at z=6 km 

4 km 

at z=6 km 

Peak Vertical 
Velocity near Flight 
Level 

Max 
10 
in s' 1 

Min 

-6 

m s' 1 

Max 

11 

m s' 1 

Min 

-4 

m s' 1 

Horizontal Scale of 
Turbulence Patch at 
Flight Level 

3 km 

3 km 

Peak <j\ n at Flight 
Level 

0.26 g* 

0.22 g 


* from Hazard Estimation Model 


Figure 4. In situ An (g’s) and o An (g’s) along path for 
Event 235_05. 



Figure 5. Three-dimensional perspective of TASS 
simulated cloud for Event 235 _05. 
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RADAR REFLECTIVITY AT 7.4KM AGL 
TASS SIMULATION OF 235-05 



Xf KM ) 


Figure 6. TASS simulated radar reflectivity (dBZ) for 
Event 235_05. 
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Figure 7. Charleston, SC NEXRAD imagery observed 
near flight path for Event 235_05. 



Figure 8. Huntsville, AL composite radar reflectivity for 
Event 232 _10. Color table for reflectivity levels are 
same as in Figure 2. 


R232 10 

Observation 

The second case occurred on 30 April 2002. In this 
event, the NASA B757 was vectored to central Alabama 
to investigate turbulence regions associated with a large 
thunderstorm complex. 5 6 Composite radar reflectivity 
from the Huntsville, AL, NEXRAD shows a line of 
convective cells that extended east to west across 
Alabama (Figure 8). Storm tops were between 10 and 12 
km (35,000 and 40,000/0 and cell motion was towards 
the east-southeast at 40 m s' 1 . Ambient winds at flight 
level, -10.5 km AGL, were from the west-northwest at 50 
m s' 1 (100 kts ), with significant vertical shear near the 
upper portions of the storm. An airmet for moderate 
turbulence had been issued indicating wide-spread 
turbulence in association with the vertical wind shear at 
the jet stream levels. The prevailing operational situation 
was an environment that precluded ‘see and avoid’ 
options for the pilots, which was the case for much of the 
region. At the flight level, low levels of radar reflectivity 
were present, often presenting no reflectivity signatures 
on the conventional ships radar-display. The 
characteristics of this event and the general operational 
hazards due to CIT are signified by the encountering of 
severe turbulence by two commercial air carriers in the 
vicinity of Event 232_10 (Figure 9). 

Numerical Simulation 

The domain size for this case is 20 km x 30 km x 
14.75 km , with periodic conditions assumed on the 
eastern and western edges of the grid. The number of 
grid points in the domain is 254 x 378 x 178, with a 
horizontal grid spacing of 80 m. The vertical grid size 
stretches with height. The domain is rotated 12 degrees 
clockwise. The initial environmental conditions were 
extracted from the Mesoscale Analysis and Prediction 
System (MAPS) forecast near time and location of event. 
MAPS is the experimental version of the Rapid Update 
Cycle. 8 

Comparisons between the observed and simulated 
environments for Event 232-10 are presented in Table 2. 
The simulated radar reflectivity field is shown in Figure 
10 and compares well with the observed radar reflectivity 
in Figure 9. 

The TASS simulation indicates overshooting tops 
with strong horizontal gradients of vertical velocity 
(Figure 11). Associated with these overshooting tops are 
regions of very weak radar reflectivity (Figure 10). 
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Figure 9. Flight path for Event 232_10. Color tables 
for reflectivity levels and aircraft RMS normal loads 
( Oj n ) along path are same as in Figure 2. Also denoted 
are two severe turbulence pilot reports that occurred 
within 30 minutes of the event. 

RADAR REFLECTIVITY AT 10 KM AGL 
TASS SIMULATION OF 232-10 CASE 
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Figure 10. Same as in Figure 6, but for Event 232 _10. 

VERTICAL VELOCITY AT 1 0 KM AGL 
TASS SIMULATION OF 232-10 CASE 


Table 2: Comparison of simulated and observed 

characteristics for Event 232_10. 


Variable 

TASS 

Observed 

Peak Storm Tops 

11 km 

10-12 km 

Peak Radar 
Reflectivity near 
Flight Level 

17.5 dBZ 

18 dBZ 
from aircraft 

Cell Motion 
(toward) 

ESE at 
35 m s' 1 

ESE at 
40 m s' 1 

Width of Convective 
Line near Ground 
(based on 20 dBZ) 

11 km 

15-25 km 

Peak Vertical 
Velocity near Flight 
Level 

Max 

11 

m s' 1 

Min 

-11 

m s' 1 

Max 

15 

m s' 1 

Min 

-21 

m s' 1 

Horizontal Scale of 
Turbulence Patch 
near Flight Level 

3-5 km 

8 km 

Peak a An near Flight 
Level 

0.38 g* 

0.45 g 


* from Hazard Estimation Model 


Wilmington , Delaware FOQA case 

Observation 

The last case is a FOQA event in which a 
commercial B-737 encountered a line of shallow 
convection on descent. An analysis of this event was 
documented by MIT Lincoln Laboratory. 9 

The airliner encountered severe turbulence on 
descent as it was vectored by ATC along a convective 
line. This event was a relatively low altitude case 
occurring at an altitude of 2.3 km. The turbulence 
encounter occurred as the airliner penetrated the leading 
edge (southeastern side) of a line of relatively shallow 
convective precipitation. 

Numerical Simulation 

The initial environmental conditions for this 
simulation were extracted from a Mesoscale Atmospheric 
Simulation System (MASS) 10 forecast near the time and 
location of the event. The domain size for this simulation 
is 14 km x 7.5 km x 6.25 km , with periodic conditions 
chosen on the lateral boundaries that are orthogonal to 
the convective line. The number of grid points in the 
domain is 283 x 154 x 128, with a uniform grid size of 
50 m. The domain is rotated 245 degrees clockwise. 



X (KM) 

Figure 11. TASS simulated vertical velocity (m s' 1 ) for 
Event 232 _10. 
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Table 3: Comparison of simulated and observed 

characteristics for Delaware FOQA case. 


Variable 

TASS 

Observed 

Peak Storm Tops 

5.25 km 

5. 5-7. 6 km 

Peak Radar 
Reflectivity near 
Flight Level 

46 dBZ 

40 dBZ 
from NEXRAD 

Cell Motion (toward) 

SE 

at 25 m s' 1 

SE 

at 20 m s' 1 

Width of Convective 
Line near Ground 
(based on 20 dBZ) 

9 km 

10-18 km 

Peak Vertical 
Velocity near Flight 
Level 

Max 
12 
m s' 1 

Min 

-10 

m s' 1 

Max 

14 

m s' 1 

Min 

-14 

m s' 1 

Peak G An near Flight 
Level 

0.33 g* 

0.36 g 

Propagation speed of 
convective line 

13 m s' 1 

6-19 m s' 1 


* from Hazard Estimation Model 


z 




Figure 12. Same as in Figure 5 but for the FOQA case. 

Comparisons between the simulated and observed 
environment are shown in Table 3. A three-dimensional 
perspective of the convective line is shown in Figure 12. 
The simulated radar reflectivity field at flight level for the 
FOQA case is presented in Figure 13. Peak radar 
reflectivity barely exceeds 40 dBZ as was true from the 
observed NEXRAD images presented by Ferris. 9 The 
associated vertical motion field (Figure 14)exhibits a 
pattern of significant updrafts and downdrafts along the 
convective line. The strongest vertical velocities are 
located on the forward side of the line and coincide with 
areas of relatively low radar reflectivity (< 30 dBZ). 



Figure 13. Same as in Figure 6, but for FOQA case. 


VERTICAL VELOCITY AT 2300 M AGL 
TASS SIMULATION OF FOQA DELAWARE CASE 



X fkm} 


Figure 14. Same as in Figure 11 but for FOQA case. 


Hazard Analyses 

The turbulence hazard is quantified for the previous 
case studies by applying hazard estimation ‘tools.’ 


Hazard Estimation Model 

As in Proctor et al. 2 , we apply hazard metric 
algorithms to the simulated data sets. 

For a particular aircraft, a An can be estimated from 
g w using aircraft look-up tables. 11 The o w fields can be 
computed for any horizontal plane in the each data set, by 
using a moving average as: 

<y w (x, y) = 

Lt 


x+— Jf+— 
2 2 


L L 

x y 


J J{ w(x',y')-w(x,y)} 2 dx'dy ' 
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for z > 9.1 km 


where the averaging interval along the x and y 
coordinates is L x , L y , respectively. The average vertical 

wind, w , is computed from the vertical wind, w, as: 


w(x, y ) 


y— | J w(x' ,y')dx'dy' (2) 

L * L > x -± 

2 y 2 


The value for the averaging interval, L x =L y = 1000 
m, is chosen to correspond to ~5 s averaging period for a 
commercial aircraft at cruise speeds. Hence, the second 
moment of the w-field is computed assuming a 1 x 1 km 
moving box. At slower speeds (such as during descent) a 
smaller averaging interval for L x and L y is assumed that 
corresponds to the shorter distance traveled by the 
aircraft in 5 s. Using formulas developed by Bowles, the 
RMS normal load can be estimated from a w for NASA’s 
B-757 by: 


y) = & w (x> y)[a(z) - b(z)log(l)] 
( 180 V 




V 


Vn(z) 


*K 


( 3 ) 


where V is aircraft’s true airspeed, W a is the aircraft 
weight in kilo-pounds, K = 1.09, and the turbulence 
length scale Z = 500 m. The parameters a(z), b(z), and 
V„( z) are: 


a(z) — C; +C2Z + C3Z +C4Z + CjZ + C(,Z + CjZ 


where: 

Ci — 0.136700 s-1 
c 2 = -4.759443e-05 mV 
c 3 - 3.453818e-08 m'V 1 
c 4 = -1.031272e-ll mV 
C 5 = 1.457209e-15 m'V 1 
Cg = -9.709422e-20 m^s 1 
c 7 - 2.448279e-24 m' 6 s 4 

b(z) = Ci +c 2 z + CjZ 2 +c 4 z 3 + c 5 z 4 + c<sz 5 


where: 

c 1 = -.0179545 s 1 
c 2 = 3.9960830e-05 m''s‘' 
Cj = - 1. 24327 le-08 mV 
c 4 = 1.760674e-12mV 
cj = -1.153713e-16mV 
c 6 - 2.807013e-21 m^s ' 1 


V n (z) = 232.21 m s' 

VJz) = Ci +c 2 z + CjZ for 3.0 km<. z < 9.1 km 
where: 

C! = 150.1493 ms 1 
c 2 = 6.2023e-3 s" ' 

Cj - 3.0298e-7 m 1 s' 1 


The accuracy of Eq. (3) has been verified from NASA B- 
757 in situ flight data. 

0 th -Order Aircraft Dynamic Model 


A second relationship for estimating a An is proposed 
that is much simpler to apply than higher-order, time- 
domain aircraft flight dynamics models. The proposed 
relationship is applied to a particular flight path of choice 
estimating the loads based on the vertical winds and is 
designed to provide an aircraft load estimate as a check 
against other techniques. The relationship has been 
verified with NASA B757 in situ data and is applied to 
the simulation results presented above. 

First, An is calculated from the horizontal gradient of 
vertical velocity along a specified flight path; i.e., : 


An = 


^—k 

8 dl 


( 4 ) 


where the running mean of w is calculated along a flight 
path / as: 


— If 2 

w(l) = — \w{V)dV ( 5 ) 

L i ik 

2 

V a is the aircraft’s airspeed, the averaging interval, Lj = 
300 m, and k is a constant determined by aircraft type and 
weight. The estimate of o An is then calculated from An 
assuming a 5s (~1 km) running window. 

Optimization of the 0 th -order model is still being 
investigated, but initial results support the constant, k = 
0.69 for the NASA B757. 

Prior to comparing the two hazard models to the 
three simulated cases, estimated loads from the 0 th -order 
model are compared to observed loads from NASAB757 
in situ data. Figure 15 and Figure 16 show comparisons 
for each of the two flight events, 235_05 and 232_10, 
respectively. The actual events compare well with the 
loads predicted by the 0 th -order model. 

For application in the simulated wind fields, the 0 th - 
order model is calculated for flight paths through each 
grid point that are parallel to either the X- or F-coordinate 
axes. Flight paths for 235_05 are calculated along the X- 
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coordinate, while flight paths for 232_10 and the FOQA 
case are calculated along the 7-coordinate. These 
headings correspond with the actual aircraft headings 
during each event. 

Figure 17-Figure 22 present the hazard analyses for 
all three events. The distributions of estimated RMS 
normal loads compare well between each of the two 
hazard models. Similarly, the estimated loads compare 
well with the measured loads from in situ (Table 4). 

Analyses of the numerical data sets show that peak 
loads from the hazard estimation models (Figure 17- 
Figure 22) are not necessarily correlated with regions of 
peak radar reflectivity (Figure 6, Figure 10, and Figure 
13). Furthermore, high turbulence intensities can occur 
outside of regions of significant radar reflectivity. 

Event 235 05 



0.00 


5.00 


10.00 15.00 

Distance (km) 


20.00 


Figure 15. Estimated and observed o An Event 235_05. 

Event 232 10 



CJ W FROM AT 7.4KM AGL 
TASS SIMULATION OF 235-05 


jr 


3 


X(KM ) 

Figure 17. Horizontal cross section of 
from the Hazard Estimation Model for 
Event 235 05. 


o An field (g’s) 
the simulated 


, FROM 0 th ORDER MODEL AT 7.4KM AGL 
TASS SIMULATION OF 235-05 


25.00 



X(KM ) 


Figure 18. Horizontal cross section of o An field (g’s) 
from the 0 th -order aircraft performance model for 
simulated Event 235_05. 

FROM ct w A T 10 KM AGL 
TASS SIMULATION OF 232-10 CASE 


0.00 10.00 20.00 30.00 40.00 

Distance (km) 


50.00 


60.00 


Figure 16. Same as in Figure 15, but for Event 232-10. 
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Figure 19. Same as Figure 17, but for Event 232_10. 
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a^FROM 0 m ORDER MODEL AT 10 KM AGL 
TASS SIMULATION OF 232-10 CASE 
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Figure 22. Same as in Figure 18, but for FOQA case. 


Summary 


Figure 20. Same as Figure 18, but for Event 232_10. 


a iB FROM cj w AT 2.3 KM AGL 
TASS SIMULATION OF FOQA DELAWARE CASE 



X (KM) 


Figure 21. Same as in Figure 17, but for FOQA case. 


Three cases of aircraft encounters with convectively 
induced turbulence have been presented. Large eddy 
simulations were conducted for each of the convective 
systems responsible for the turbulence and their results 
compare well with observations. Two hazard metric 
models have been introduced and validated with aircraft 
in situ data. The models’ predicted loads compare 
reasonably well with the observed aircraft loads. The 
LES data sets and hazard metric models can be useful 
tools for characterizing CIT turbulence and evaluating 
turbulence detection sensors. 
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Table 4. Comparisons between observed and modeled o An for each of the three cases. 


Source 

Peak c An (g’s) 

Event 235_05 

Event 232_10 

FOQA Case 

In situ 

0.22 

0.45 

0.37 

Hazard Estimation Model (g w ) 

0.26 

0.38 

0.33 

0 th -Order Aircraft Perf. Model 

0.27 

0.43 

0.45 
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